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Following the successful application of the relativistic heavy quark(RHQ) action to the charmo- 
nium and charm-strange meson spectrum, we study here the bottom system to explore the validity 
of this method in a regime with larger heavy quark momenta. The spectrum is calculated using 
the same dynamical 2+1 flavor, 24 3 x 64 domain wall fermion lattice configurations generated by 
the RBC and UKQCD collaborations and used in the earlier charm study. The 3 parameters in the 
RHQ action are determined by matching to the experimental bottom-strange meson masses and 
extrapolated to the chiral limit from three different sea quark masses. We predict the bottomo- 
nium mass spectrum and compare it with experiment. A theoretical estimation is also carried out 
to understand the 0(a 2 p 2 ) systematic errors found in the numerical study. 
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1. Introduction 

Heavy quark physics is a long standing challenge for lattice QCD due to the fact that the 
masses are large in current lattice units. Here ma <C 1 is no longer true and terms containing (ma) n 
(with a the lattice spacing) become significant. This problem must be addressed by using effective 
field theories since a direct simulation is beyond the current computation power. Various heavy 
quark effective actions were developed for different physical systems. See Refs. [jl|, ||, BL |4j] for 
reviews on this topic. 

In Ref. []|] we studied the relativistic heavy quark(RHQ) action ^ 0, ||, ^] on charmonium 
and charm-strange systems, and found excellent accuracy for this application on charm quarks. We 
determined the RHQ parameters at the 1% level, predicted XcO and % c \ masses with sub-percent 
accuracy, and determined the lattice scale on the j8 = 2.13 , 24 3 x 64 lattice ensemble with an 
accuracy at least as good as that of other methods. In this proceeding we will continue to explore the 
three-parameter RHQ method in bottom systems, a regime with much larger heavy quark momenta. 
We start with the bottom-strange system to determine the RHQ parameters as it has smaller a 2 p 2 
discretization errors since p ~ Aqcd, and then calculate the bottomonium states. To our surprise 
the calculations on bottomonium states have errors around 30 MeV compared to the experiment. 
An attempt to theoretically estimate the 0(a 2 p 2 ) errors expected in this numerical study is also 
carried out at the end. 

We start by briefly introducing the framework we have been using before. The lattice form of 
the action, following the formulation proposed in M, m, can be written as: 



5 = m + %D -^Oo + C f-D-^a(D) 2 - a£ ^c P cr MV iy v \ V n , (1.1) 



In heavy quark system, the temporal covariant derivative Do is on the order of ma but the spatial 
derivatives D, are of order Aqcd® or a s ma depending on the system under investigation. We found 
that only three free parameters are needed in the action: mo, cp and £. If the parameters are 
correctly tuned, the action will have small cutoff effects: 0(Kqcdo) 2 for heavy-light systems and 
0(a s ma) 2 for heavy quarkonium. 

The lattices used in this work are the dynamical 2+1 flavor 24 3 x 64 DWF lattice configurations 
generated by the RBC-UKQCD collaborations. For better statistics, we place sources at different 
time slices separately for each configuration. Binning the data does not give any significant change 
in the size of the errors, indicating the auto-correlation of the lattice configurations is negligible. 



volume 


L s 


(m sea ,m s ) 


Traj(step) 


# of configs 


sources 


24 3 x 64 


16 


(0.005,0.04) 


900-6880(20) 


300 


0,32 


24 3 x 64 


16 


(0.01,0.04) 


1460-5060(40) 


91 


0,16,32,48 


24 3 x 64 


16 


(0.02,0.04) 


1885-3605(20) 


87 


0,16,32,48 



2. Determination of the RHQ action 

To determine the action in a way in which errors are controlled, the parameters are fixed 
by matching physical observables, sensitive to those parameters, to their experimental or theo- 
retical values. These parameters are determined for each ensemble with three different light sea 
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quark masses and extrapolated to the chiral limit. In addition to mass combinations of pseudo- 
scalar (rjb,B s ), vector (T,B*), scalar (%w) an d axial- vector (Xbi) mesons in heavy-heavy (hh) and 



heavy-strange (hs) systems (10p, we also calculate the tensor state (hb) for heavy-heavy system. 
Specifically we compute: 

• Spin-averaged: mf a = \(m nb + 3m T ), mf a = \{m Bs + 3m B *) 



Hyperfme splitting: m n h n s = m r - m nb , m n h l s = m B * - m Bs 
Mass ratio: where E 2 = m\ + ^p 2 , mi: rest mass, 1112: kinetic mass. 



1112 



• Spin-orbit averaged and splitting (hh): m^ s = m Xbj 



- m Xb0 , m^g a 



i( m Xb + 3m Z 



Heavy-heavy l Pi state hb 



Here are some sample plots of the bottomonium^) and bottom-strange(B s ) pseudo-scalar corre- 
lators. One should notice the correlators are falling about 70 orders of magnitude for bottomonium 
and about 40 for bottom strange. We must be careful that the heavy propagators are sufficiently 
accurate when the time is large. 
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Figure 1: sample correlators of bottomonium and bottom-strange pseudo-scalar. 

The matching method is the same as used in Ref. We again use the linear ansatz relat- 
ing the three parameters (Xrhq) and the corresponding measured quantities (Y(a)) and determine 
the parameters by minimizing the % 2 defined in that reference. The linear approximation only 
holds in a certain region of the parameter space, which we estimate from our earlier work using a 
16 3 x 32 lattice. The linear coefficients are calculated directly using finite differences from seven 
point Cartesian parameter sets, which are centered at {7.3,4.0,4.3} and with extent {0.5,1.0,0.3}. 
We examine the linearity of the parameter dependence in this region and required that the finally 
matched RHQ parameters actually lie within the region examined. 



2.1 Numerical details 

For all the calculations we use a box source with size 4, which is not optimized to give the best 
plateau. The physical strange quark mass used in the bottom-strange spectrum is 0.034 in lattice 
units [p|]. Some sample effective mass plots for different states are shown for the m sea = 0.005 
ensemble, Fig |2| and ||(left). We can see from the plots that the plateau for XbO>%b\ an d hb states 
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are not quite clear so their masses are subject to more systematic errors. All meson correlators are 
fit to a single state with exponential decay. The fitting ranges are chosen from an examination of 
the plateau in the effective mass plots, and they are summarized below in Tab. [T]. 



fitting range 


ribTC 


XwlXbilhb 


BJB* 


m sea = 0.005 


[14,30] 


[5,12] 


[10,25] 


tn S ea = 0.01 


[11,30] 


[5,9] 


[10,25] 


m sea = 0.02 


[11,30] 


[6,15] 


[13,18] 



Table 1: The time fitting ranges for different mesons at different sea quark masses. 

The momentum dependence is studied for both the 7]b and T mesons and the mass ratio m\/m,2 
extracted from each are quite consistent. We use results from the T momentum dependence with the 
three lowest momenta. Please note that although we are using bottom-strange spectrum to fix the 
RHQ parameters, the mass ratio m\/m2 is determined from bottomonium because it is much better 
determined than it is from bottom-strange. To ensure the precision of the heavy quark propagator, 
we use an extremely tight CG stopping condition of 10 -60 . The CG stopping conditions are set to 
10~ 8 for light quark propagators as before. 
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Figure 2: effective masses for the mesons on m sea — 0.005 ensembles: (left)r//, and T; (right )B S and B* 
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Figure 3: (left) Effective masses for the mesons: XbOtXbi an d hb- (right) the prediction of rjt, mass in the 
chiral limit compared with the experiment. 
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3. Analysis and Results 

Again we list all the quantities calculated in this work: (1)^ (2)ms* — (3)i(mB s + 3ms*) 



(4)mx 



m Vt ( 5 )? K?* + 3m r) {6)m XM - m XbQ K , 



(7)1 (m, 



">2 



3m y J (8)m/, 



3.1 Heavy-strange sector: the RHQ action parameters 

As we mentioned before, except for the mass ratio m\/m.2 which is determined from the T 
momentum dependence, only the heavy-light, B s and B* masses are used to fix the RHQ action 
parameters for the bottom system. Throughout the analysis, the lattice scale is taken to be a -1 = 
1.73 GeV [O]. From matching to quantities (l)-(3), the RHQ parameters and the corresponding 







cp 


C 


0.005 


7.37(7) 


3.84(40) 


4.21(3) 


0.01 


7.28(9) 


3.28(40) 


4.21(3) 


0.02 


7.30(11) 


3.52(53) 


4.24(4) 


-m res 


7.38(12) 


3.93(54) 


4.19(4) 



Table 2: The RHQ action parameters determined from matching the quantities (l)-(3), and exttapolated to 
the chiral limit. 



3.2 Heavy-heavy sector: the predictions on bottomonium states 

What would we expect for the size of the errors on the bottomonium masses computed in 
this way? A naive estimate would be the size of a typical order a 2 operators, for example, O = 
*Py D D 2x ¥. The typical velocity in bottomonium system v ~ 0.1, which is determined by the 
splitting T(2S) -T(1S) ~ mv 2 ~ 500 MeV. This indicates (6) ~ ^ ~ 300MeV. 

Let's turn to the numerical results given by this work. All calculations are made from a com- 
plete calculation of each jackknife block from which the average and errors are determined. The 
results for various quantities are summarized in Tab. 3.2 and compared to experimental values and 
other lattice calculations. Our results include only statistical errors. Experimental results are from 
PDG unless otherwise specified. We can see that most of our calculated bottomonium masses are 
within 30 MeV(<l%) and one standard deviation of the experimental value. Of course, the results 
for the mass splittings give a more precise test: the spin-orbit splitting from our calculation is about 
2.5 standard deviations from the experiment and the hyperfine splitting is further away from the 
value newly determined by BaBar collaboration [13] or lattice NRQCD calculation [14]. But the 
accuracy of our results indeed comes as a surprise! We were expecting errors of hundreds MeV but 
only see errors several times smaller. 

A possible better theoretical estimate would use a simple hydrogen-like Coulomb model to 
describe bottomonium. In that case, we could define a T state as: 

lT ' mj) = I d -^0^I<^\ p -^ 

where a and b are free field quark and anti-quark annihilation operators as in 

m = f 1 l g {u s (p)e i ^a{p, S )+v s {p)e- i ^(p, S )} (3.1) 
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quantities 


RHQ(MeV) 


Exp. (Me V) 


NRQCD(MeV) 


m Vb 
ntx 

m Xw 
m XM 
m h 
ntx ~ m nb 

m Xbi ~ m Xbo 


9420(14) 
9444(17) 
9873(15) 
9897(16) 
9908(17) 
23.7(3.7) 
24.0(3.5) 


9389(3)(3) [13] 
9460 
9859 
9893 

71(3)(3) [13] 
33.34 


9900(6) [14] 
61(14) [14] 



Table 3: The predictions on individual masses and splittings, compared to experiment and lattice NRQCD 
calculations. Results from our calculation include statistical errors only. 



and <p(p) is the Fourier transform of the hydrogen atom wave function with mass m/2, the reduced 
mass of the bottomonium system. 

The size of the matrix element of the operator O can then be estimated and yields the result: 



(O) 



J 3 4 2 
~ o m b a s a 



40MeV m b = 4.0GeV, a s = 0.25 
U6MeV m b = 4.0GeV, CC S = 0.35 



Thus the size is very sensitive to the strong coupling constant a s . As an accurate a s is not available 
we can not make any strong statement from this approach. 
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Figure 4: The mass predictions in the chiral limit compared with the experiment, (left) T; (right) 



4. Conclusion 

We have studied the bottomonium and bottom-strange system. We calculate the RHQ pa- 
rameters by matching the bottom-strange system to the experiment at each sea quark mass and 
extrapolate to the chiral limit. The bottomonium spectrum is then determined in the chiral limit 
within 30 MeV of the experimental values. A naive estimate of the size of a a 2 p 2 operator sug- 
gests the error expected in bottomonium spectrum should be a few hundred MeV, while a more 
careful hydrogen-like Coulomb model suggests the error size is very sensitive to the strong cou- 
pling constant a s . Other phenomenological models might be useful to estimate the error size more 
accurately. In summary, our application of the RHQ action on the bottom system yields surpris- 
ingly accurate results. Our possible next steps include study of heavy-light system which could 
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Figure 5: The mass predictions in the chiral limit compared with the experiment, (left) Xbol (right) Xbi- 

double-check the validity of the RHQ method in this regime, more calculations on other states like 
be system and nucleons with one or more heavy quarks, and matrix elements. These calculations 
could also be repeated using the new 32 3 x 64 lattices to better determine the continuum limit. 
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